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Abstract. We consider the geometric phase and quantum tunneling in vicinity of 
diabolic and exceptional points. We show that the geometric phase associated with 
the degeneracy points is defined by the flux of complex magnetic monopole. In weak- 
coupling limit the leading contribution to the real part of geometric phase is given 
by the flux of the Dirac monopole plus quadrupole term, and the expansion for its 
imaginary part starts with the dipolelike field. For a two-level system governed by the 
generic non-Hermitian Hamiltonian, we derive a formula to compute the non-adiabatic 
complex geometric phase by integral over the complex Bloch sphere. We apply our 
results to to study a two-level dissipative system driven by periodic electromagnetic 
field and show that in the vicinity of the exceptional point the complex geometric phase 
behaves as step-like function. Studying tunneling process near and at exceptional 
point, we find two different regimes: coherent and incoherent. The coherent regime is 
characterized by the Rabi oscillations and one-sheeted hyperbolic monopole emerges 
in this region of the parameters. In turn with the incoherent regime the two-sheeted 
hyperbolic monopole is associated. The exceptional point is the critical point of the 
system where the topological transition occurs and both of the regimes yield the 
quadratic dependence on time. We show that the dissipation brings into existence of 
pulses in the complex geometric phase and the pulses are disappeared when dissipation 
dies out. Such a strong coupling effect of the environment is beyond of the conventional 
adiabatic treatment of the Berry phase. 
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1. Introduction 

The geometric phase was first discovered by Berry pQ in the context of the adiabatic 
cyclic evolution of the quantum system. Later Wilczek and Zee [2] generalized the 
Berry's phase allowing the transported states to be degenerated, and Aharonov and 
Anandan [3] extended Berry's result to non-cyclic and non-adiabatic variation of the 
parameters of the Hamiltonian. However, while notion of the geometric phase for pure 
states is well defined, the definition of a geometric phase in open quantum systems is 
still an unsolved problem. 

The first important step towards to the consistent description of geometric phase 
for open systems was done by Garrison and Wright [I] . They removed the restriction to 
unitary evolution considering the quantum system when the time evolution in governed 
by non-Hermitian Hamiltonian. Beginning with the classical works by Weisskopf and 
Wigner on metastable systems [3 modeling the dissipative quantum systems by 
effective non-Hermitian Hamiltonian is well known. It was observed that for the 
system being initially in the metastable state ^(O), its evolution is described by the 
effective non-Hermitian Hamiltonian as follows: ip(0) — > ip(t) = e~ tHeft ip(0)+ decay 
products. Its derivation can be made by separation of the full Hilbert space into the 
intrinsic discrete part and continua and applying the projection operator technique 
to eliminate the continuum part (see e.g. [3 El [9]). Since the Garrison and Wright 
paper has been published, the geometric phase for quantum systems governed by the 
non-Hermitian Hamiltonian and complex-valued geometric phase effects in dissipative 
systems were studied by various authors (for discussions and references see, e.g., 



Actually, for correct description of open quantum system a density-matrix approach 
is required. In particular, when the state of the system is changed negligible during 
the time scale characterizing the decay of the reservoir correlation function, the Born- 
Markov approximation is valid, and an open quantum system can be described by the 
following general master equation written in the Lindblad form {fi = 1) [25J: 



where H is the Hamiltonian of the uncoupled system and r& are transition operators 
describing the coupling to the reservoir. The commutator of the density operator p 
with a Hamiltonian H represents the coherent part of evolution and the remaining 
part corresponds to the decoherence process on account of the interaction with an 
environment. The first analysis of geometric phase based on the master equation has 
been done by Ellinais et al [26] and Gamliel and Freed [27] . 

To make connection with the non-unitary evolution of quantum system, consider 
the quantum jump approach to open systems, in which the evolution of system may be 
separated in a part describing "no-jump" and spontaneous jumps part as follows. The 
total evolution time T is divided into a sequence of intervals At = T/N. Then, the 
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jump-free trajectory being governed by the effective non-Hermitian Hamiltonian 



is interrupted by instantaneous jumps caused by jump operators T k : p{t m ) — ► p(t m+1 ) = 
E n k= Q W k p(t m )Wl where W = 1 - iHAt, W k = v^AtT k , for k = 1,2, ... ,n, and 
t m = mAt [28] [29] . Recently, these ideas has been employed to show how the geometric 
phase can be modified for open systems [2S] 13T)] . 

Recent experimental results providing evidence for the 'magnetic' monopole in 
the crystal-momentum space [31] and emerging of 'fictitious magnetic monopoles' 
in the anomalous Hall effect of ferromagnetic materials, magnetic superconductors, 
trapped A-type atoms, anisotropic spin systems, noncommutative quantum mechanics, 
in ferromagnetic spinor Bose-Einstein condensate, etc. [HH [32], [33] EU [35] [36] [37] 138] . 
has been caused a rebirth interest in the Dirac monopole problem. 

The 'fictitious' monopoles appear in the context of the Berry phase as follows. 
Assume that for adiabatic driving quantum system the energy levels may cross [1]. 
Then, in the commonest case of double degeneracy with two linearly independent 
eigenvectors, the energy surfaces form the sheets of a double cone, and its apex is 
called a "diabolic point" [39]. Since for generic Hermitian Hamiltonian the codimension 
of the diabolic point is three, it can be characterized by three parameters R = 
(X,Y,Z). The eigenstates \n, R) give rise to the Berry's connection defined by 
A„(R) = i(n, R| Vr|ti, R), and the curvature B n = Vr x A n associated with A n is 
the field strength of 'magnetic' monopole located at the diabolic point [U [10] • The 
Berry phase 7„ = § c A n ■ dH is interpreted as a holonomy associated with the parallel 
transport along a circuit C [40J. 

In opposition to the diabolic point, at which the eigenvalues coincide while the 
eigenvectors still remain distinct, an exceptional point occurs when the eigenvalues 
and eigenvectors coalesce. This type of degeneracy is associated with open quantum 
systems and non-Hermitian physics [28] [29] UH [TT] H2] . The exceptional points have 
been observed in various physical systems: laser induced ionization of atoms, microwave 
cavities, "crystals of light", in optics of absorptive media, electronic circuits, etc. 



In this paper we consider the geometric phase and tunneling process near and 
at diabolic and exceptional points. We show that for general non-Hermitian system 
the geometric phase associated with the degeneracy is described by complex magnetic 
monopole. We find that the exceptional point is the bifurcation point of the complex 
geometric phase in the parameter space and the real part of the latter has a jump 
discontinuity at the exceptional point. We show that the exceptional point is the critical 
point of the quantum-mechanical system, where the topological phase transition in the 
parameter space occurs. 

Studying the tunneling process in the vicinity of exceptional point we found two 
distinct regimes: coherent and incoherent. The coherent tunneling is characterized by 
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the Rabbi oscillations, also known as quantum echoes. We also show that the dissipation 
brings into existence of pulses in the real part of the geometric phase. Such a strong 
coupling effect of the environment disappears in the absence of dissipation. 

2. General results on behavior of the eigenvectors at diabolic and 
exceptional points 

As known, in parameter space, a set of exceptional points defines a smooth surface of 
codimension two for symmetric/nonsymmetric complex matrix, codimension one for a 
real nonsymmetric matrix, and exceptional point's do not exist for real symmetric or 
Hermitian matrix [IB]- Let H(X) be a complex N x N matrix smoothly dependent on 
m real parameters Aj, i runs from 1 to m. For Afc(X) being the eigenvalues of H(X), 
we denote by \^k(X)) and (ipk(X)\ be the corresponding right/left eigenvectors: 

H\tpk) = A fc |V>*), (i>k\H = X k {i>k I (3) 
Both systems of left and right eigenvectors form a bi-orthogonal basis 



y\^0tA = 1 (^|^) = , k^k>. (4) 

Using the decomposition of unity (j4]), we obtain 

|V) = J>^>, (^l = $>^l (5) 

i i 

where 

a, = — , /5j = — . (b) 

We assume that exceptional point occurs for some value of parameters X = X c . At 
the exceptional point the eigenvalues, say n and n + 1, coalesce: A„(X C ) = A n+ i(A" c ), 
and the corresponding eigenvectors coincide, up to a complex phase, yielding a single 
eigenvector \i/jep)- Now, applying (j4]) for k = n and k = n + 1 we find that at the 
exceptional point the normalization condition is violated, (iPepI^ep) — 0. This leads to 
the serious consequences for the global behavior of the states on parameter space. 

Since at exceptional point both eigenvalues and eigenvectors merge forming a 
Jordan block, it is convenient to introduce the orthonormal basis of the related invariant 
2-dimensional subspace as follows: 

(n\n) = l, (n + l\n + 1) = 1, (n|n + l)=0 (7) 

Assuming that all other eigenstates are non-degenerate, we find that the set of vectors 
{\Xk), (Xk\,} , where \\ n ) = \n), (Xn+i\ = (n + 1|, and 

|») = -J#=, (U = -^=, for k * n, n + 1 

form the bi-orthonormal basis. Using this basis we expand an arbitrary vector ip as 

\^) =Y, c ^ X )\Xk{X)) (8) 
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with the coefficients of expansion defined as = (xk\4>)- 

From the orthogonality condition, one can see that if \i/j(X)) — > \ipEp) while 
X — > X c , all coefficients Ck {k ^ n,n + 1) vanish at exceptional point. Thus in the 
neighborhood of exceptional point only the terms related to the invariant subspace 
make substantial contributions and the iV-dimensional problem becomes effectively two- 
dimensional [HI [50]. Similar conclusion is valid for the diabolic point, excepting that 
at the diabolic point the eigenvectors remain distinct. 

Now, let being the eigenvector of H, 

Hmx)) = \(x)mx)), (9) 

belongs two-dimensional invariant subspace defined by the Jordan block. This implies 
that in the expansion (jSJ) for k ^ n, n + 1 all Ck = 0, and we may write 

\i/;)=a\n)+P\n + l) (10) 

For X ^ X c , one has two eigenvalues A + and A_ and the corresponding eigenvectors 
\if)±) are given by 

\ip ± ) = a±\n) +(3±\n+ 1) (11) 

The coefficients A±,a± and (3± are found from the two-dimensional eigenvalue 
problem 





(12) 



an ai2 

where the matrix elements are determined as 

On = (n\H\n), a\2 = (n\H\n+ 1), 
a-21 = (n+ l\H\n), a 22 = (n + l\H\n + 1) 

It is convenient to introduce the following notations: 

a u + a 2 2 v a u + a 21 , , 

A o = ^ ' ~ 2 ' 

y _ a 21 — Ol2 „ _ CLu — 2 2 (lA) 

2i ' ~ 2 ' [ ' 

Then reads 



\ + Z X- 
X + iY A 





(15) 



Solving the characteristic equation for ffl5l) . we obtain 

A± = A ± VX 2 + Y 2 + Z 2 (16) 

The eigenvalues coalescence at the point R = 0, that yields the diabolic point if 
X = Y = Z = 0, and the exceptional point otherwise. The detailed study of two- 
dimensional problem will be presented in the following sections. 
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3. Degeneracy, geometric phases and complex 'magnetic' monopoles 

Following [I], let us consider the time dependent Schrodinger equation and its adjoint 
equation: 

i^Mt)) = H(X(t))Mt)), (17) 

^(9{t)\ = m)\H(X{t)), (18) 

where H is the non-Hermitian Hamiltonian. 

Let (ip n (X)\ and \ijj n (X)) be left (right) eigenstates corresponding to the eigenvalue 
E n , then in adiabatic approximation the complex geometric phase is given by [H fl2lfT3] 

= . / mmux)) 

L (MX)\MX)) 

generalizing Berry's result to the dissipative case. Further we assume that the 
instantaneous eigenvectors form the bi-orthonormal basis, (ip m \ip n ) — 8 mn . This can 
alter the geometric phase f|T9|) up to the topological contribution irn, n G Z [511 [52] . 

For a non-Hermitian Hamiltonian, validity of the adiabatic approximation is defined 
by the following condition: 

^ | $ m \dH/dt\MX)) ^ i (0C) ^ 

This restriction is violated nearby the degeneracies related to any of diabolic point or 
exceptional point, where the eigenvalues coalesce. 

Since the adiabatic approach cannot be applied in the neighborhood of degeneracy, 
we will consider non-adiabatic generalization of Berry's phase introduced by Aharonov 
and Anandan [3] and extended by Garrison and Wright to the non-Hermitian systems 
as follows [4J. Let an adjoint pair {|*(t)}, being a solution of Eqs. (fT7jh ffTgD 

satisfies the following condition 

|*(r))=exp(«V)|*(0)> J (21) 

(*(T)|=expH»<$(0)|, (22) 

where ip is complex, and a modified adjoint pair such that 

\ X (t)) = eMif(t))m)), (23) 

where f(t) is any function satisfying f(t+T) — f(t) = <p(t). The total phase ip calculated 
for the time interval (0, T) may be written as (p = 7 + 5, where the "dynamical phase" 
is given by 

S=~ f \x{t)\H\x{t))dt. (25) 







and for the geometric phase 7 one has. 



7 = * f \x{t)\^X{t))dt (26) 
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This yields the connection one-form and the curvature two-form as follows 

A = i(x\d X ), F = dA. (27) 

Note that real part of geometric phase (l26i) besides of the usual Berry phase, contains the 
contribution of environment. Its imaginary part changes the amplitude of the density 
matrix and implies mixture of the initially pure states. 

Geometric phase 7 for an arbitrary quantum evolution can be obtained, also, from 
the total phase 74 by subtracting the dynamical phase 7c z [23] : 

7 = It- Id, (28) 

where 

7t = arg<¥(0)|*(*)), and ld = -i J (*(t)\^\*(t))dt. (29) 

We adopt and generalize this definition of the geometric phase for non-Hermitian 
quantum evolution as follows (see also [5T]): 

,4,„|ai + ,f, i|t) A lffl4 PC, 



2 (*(0)|*(t)) Jo dt> 
As can be observed, (j3"01) yields gauge invariant definition of the geometric phase with 
respect of gauge transformations: 

|#) e te |*), (*| -> e _te (*|, a G C. (31) 
3.1. Two -lev el system and 'magnetic' monopoles 

As has been mentioned before, in vicinity of degeneracy point behavior of iV-dimensional 
system can be described by the effective two-dimensional quantum system. In what 
follows, we consider in detail the complex geometric phase associated with the generic 
non-Hermitian Hamiltonian: 

*=(£ + + £ w)' wec (32) 

X, Y, Z G C being complex parameters. 

For the Hamiltonian ([3*21 the exceptional point is determined by equation 

X 2 + Y 2 + Z 2 = 0, (33) 

and defines a hypersurface of complex codimension 1 in C 3 , which also can considered as 
a smooth surface of codimension 2 in 6-dimensional real space. Note, that the diabolic 
point being just a point in 3-dimensional complex space C 3 , is located at the origin of 
coordinates. 

The solution of the eigenvalue problem 

H\u) = \\u), (u\H = \(u\, (34) 

where \u) and (u\ are the right and left eigenvectors, respectively, is given by 

A± = A ± R, (35) 
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where 



/ fl+Z . 9 R-Z 

COS 2 = V^' Sm 2 = V^' (38) 

Vi? 2 - Z 2 VR 2 -Z 2 y ' 

and 6, (p are the complex angles of the complex spherical coordinates: 

X = RsinO cosip, Y = Rsm8sin(p, Z = Rcos6 (40) 
Finally, for R ^ the following relationships hold 

(w±|u T )=0, (u±|u±) = l (41) 

As seen from Eq. (I35j) . the coupling of eigenvalues A + and A_ occurs when 
X 2 + Y 2 + Z 2 = 0. This yields two cases. The first one, defined by 9 = 0, — 0, 

yields two linearly independent eigenvectors. The point of coupling is known as the 
diabolic point, and we obtain 



u + )=[ n }, (u + \ = (l,0), \uJ}=[ 1 ,(w_| = (0,1), (42) 





The second case is characterized by coupling of eigenvalues and merging of eigenvectors, 
as well. The degeneracy point is known as the exceptional point, and we have: 
\u + ) = e tK \uJ) and (u + \ = e~* K (£t_|, where k G C is a complex phase. Hence, the 
violation of the normalization condition (|4l|) is occurred at the exceptional point, and 
we have (u±\u±) = 0. 

Let us assume that the exceptional point is given by Ro = (X ,Y , Z ). Then, if 
Z ^ 0, using Eqs. (EEJ) - (SO]) we obtain 

tan| = ±i, e»» = ^l±^, (43) 

and, thus, at the exceptional point — > ±00. In turns, if Z = 0, we obtain X = ±iY Q . 
This implies 8q = tt/2, and *<s<p — > ±00 at the exceptional point. 

Inserting formulae (I3"6"|) - (I3TI) for \u±) and (u±\ into (12T1) . we obtain the connection 
one-form as 

A = q(l- cos 6)d<p (44) 

where q = =Fl/2 and upper/lower sign corresponds to \u±), respectively. The related 
curvature two-form reads 

F = dA = q sin 9 d6 A dip, 6,peC, (45) 
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and in the complex Cartesian coordinates the connection one-form and the curvature 
two-form can be written as 

F=±e ijk X k dX i AdX*. (47) 

The obtained formulae describe complex "magnetic monopole" with a charge q and the 
field B = *F given by 

B = q^, (48) 

where R = (X,Y,Z), and X,Y,Z £ C. The field of the monopole can be written as 
B l = —d&/dX\ where the potential $ = q/R. 
Computation of geometric phase yields 

7 = j A, (49) 

where integration is performed over the contour C on the complex sphere S 2 . Applying 
the Stokes theorem we obtain 

1 = ^F = qQ(C), (50) 

where £ is a closed surface with the boundary C = <9£, and Q,(C) is the complex solid 
angle subtended by the contour C. 

Generally, the complex magnetic monopole is emerged in quantum mechanical 
systems admitting SL(2, C) symmetry. In particular case of 50(3) symmetry the related 
degeneracy is referred as the diabolic point, and the formula ( 1451) reproduces the classical 
Berry result on two- fold degeneracy in parameter space pp. For the exceptional point 
the field of the corresponding "monopole" represents a complicated topological charge 
rather than a point-like magnetic charge. In what follows we discuss two particular 
applications: Hyperbolic monopole and Complex Dirac monopole. 

Hyperbolic monopole. Let us consider the following non-Hermitian Hamiltonian: 

H=( A ° + fa x ~ iy ), *,y,zeu (51) 

\ x + ty \q — iz J 

1 /2 

The eigenvalues of the Hamiltonian given by X± = Xq ± R, where R = (x 2 + y 2 — z 2 ) , 
coalesce at the point R = 0. In addition, the exceptional point is represented by the 
double cone with the apex at the origin of coordinates, and the diabolic point is just 
located at the origin of coordinates. 
Applying ( 1461) - (1471) . we obtain 



q(xdy-ydx) iq 
p , p , ■ v , and F = -^e ijk ~ 
R(R + iz) R A 



A = „/r> \ ' anci F = ~^£ijkX dx l A dx J , (52) 

This yields 



B = iq R^' R= ( x '^' z ) ( 53 ) 
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Hence, the field B can be written as = — V$, where $ = —iq/R. 

For R 2 > 0, the surface defined by R = const is the one-sheeted hyperboloid, which 
can be considered also as the coset class SU(1, ~ SO(2,l)/SO(l,l) [HU 155] . 
Introducing the inner coordinates (8, <p) as 

x = R cosh 9 cosy?, ?/ = i? cosh 9 sin yj, z = .Rsinh#, (54) 

where — oo < 9 < oo, < < 27r, we obtain 

A = q(l — i sinh 8) dip F = —iq cosh 8 d8 A dip (55) 

The obtained one-sheeted hyperbolic monopole carries the imaginary total charge iq 
and, in contrast to the point-like Dirac monopole, has the singularity on the surface of 
the double cone identified with the exceptional point. 

For R 2 < 0, the surface characterized by z 2 — x 2 — y 2 = const is the two-sheeted 
hyperboloid. A convenient parametrization is given by 

x = R sinh 8 cos ip, y = R sinh 8 sin ip, z = Rcosh.8 (z > 0) (56) 
x = i? sinh # cos y = R sinh 8 sin ip, z = —Rcosh.8 (z < 0) (57) 

where < 8 < oo, < ip < 2n, R = (z 2 — x 2 — y 2 ) 1 ^ 2 , and we preserve the 
same notations for the angular coordinates, as above. The relevant coset class is 
SU{1, l)/U(l) ~ SO(2, l)/50(2), and once return to fl52]) and ((53), we obtain 

A+ = qr(l-cosh0)d^, F + = -gsinh^c/^ Arfy? (2 > 0), (58) 
A~ = g(l + cosh8)d(p, F~ = q sinh 8 d8 A dip (z < 0). (59) 

Hence, the obtained monopole carries a real total charge given by —q. Note, that 
the two-sheeted hyperbolic monopole has been already appeared in the literature in 
connection with the geometric phase (see, e.g. [HU ES, 56J). Moreover, as has been 
pointed out by Jackiw [55], this is a topologically trivial case, and the curvature may 
be removed by a globally well-defined canonical transformation. 

As can be easily shown, in the case of the one-sheeted monopole the corresponding 
potential $ is imaginary, and for the two-sheeted hyperbolic monopole $ is a real 
function. In Fig. CD the surfaces of 3$ = const and 3ft$ = const related to one-sheeted 
and two-sheeted hyperbolic monopole, respectively, are depicted. 

The hyperbolic monopoles are appeared in a wide class of physical systems 
admitting 50(2, 1) invariance (for discussion see, e.g. [HUEHIEE] and references therein). 
For instance, the hyperbolic monopole is emerged in a two-level atom interacting with 
an electromagnetic field (see Section 4.1.3) 

Complex Dirac monopole. Let us consider the non-Hermitian Hamiltonian written as 

/ A + z-ie x-iy \ 
H = \ , x,y,zeR. 60 

\ x + ty A — z + is I 

The computation of the 'magnetic' field B yields 

B =§> R=(x,y,z-ie), (61) 
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Figure 1. Hyperbolic monopole. Left panel: R 2 > 0, one-sheeted hyperbolic 
monopole. The surfaces 9<I> = const are depicted. Right panel: R 2 < 0, two-sheeted 
hyperbolic monopole. The surfaces Sft^ = const are presented. The exceptional point 
is realized as the double-cone (not plotted). 

where R = (x 2 + y 2 + z 2 — e 2 — 2iez) 1 ^ 2 . The exceptional point obtained as the solution 
of equation 

x 2 + y 2 + z 2 - e 2 + 2iez = 0, 

is the circle of the radius e a the plane z = 0. 

The field of the monopole, B = — V$, is defined by the complex potential $ = q/R, 
depicted in Fig. El Setting r = (x 2 + y 2 + z 2 ) 1 / 2 and using the real spherical coordinates 




Figure 2. Complex Dirac monopole. The surfaces 5R$ = const (left) and 3$ = const 
(right) are plotted. The exceptional point is depicted by the circle of the radius 
r = l{e = l). 

(r, a,P), we have R = [r 2 — 2iercosa — e 2 ) 1 ^ 2 . Then, we may expand the potential $ 
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as follows: 



• 2 - 2ie 

For r>e, this yields 



$ = 7 2 = =j = ?E ^(cos«) (62) 

V r 2 - 2zer cos a - e 2 ~ r l+1 



q pcosa Q(3cos 2 a — 1) , , 

® = -+*-^-^^ - + •••> (63) 

where q is the monopole charge, p = qe is the dipole moment, and Q = qe 2 is the 

quadrupole moment. 

The geometric phase of the ground state is found to be 

f ( f cos a — ie \ , „ . . 

1 = q l = - W (64) 



., , \/r 2 — 2ier cos a — £ 2 • 

Using multipole expansion (I63I) . we obtain for 7 the following expression: 

. / psin 2 a / 3Q sin 2 a cos a in 
7 = 1m + %J — — d(3 - j> 2r2 + . . . , (65) 

where jm = q § c i}- — cos a)d(3 is the contribution of the Dirac monopole at the origin, 
the second term describes the dipole contribution to the imaginary part of the geometric 
phase and the third term the quadrupole contribution to its real part. 

Let us consider the closed curve C parameterized by j3 with the complex angle 
6 = const. Then for q=l/2 the geometric phase f l64l) becomes 

1 = 7r { 1 ~^+(z-is) 2 ) (66) 

where p = (x 2 + y 2 ) 1 ^ 2 . 

A complex Dirac monopole and related complex Berry's phase appear in wide class 
of open systems, where the Hamiltonian 

h = B(t) ■ a - ^r f r 

includes spontaneous decay T = yfea- as a source of decoherence (see, e.g., |24J and 
references therein). For instance, it emerges in a two-level atom driven by periodic 
electromagnetic field E(t) = $t(£(t) exp(iut)), with £(t) being slowly varied, as follows. 
In the rotating wave approximation, after removing the explicit time dependence of the 
Hamiltonian with a suitable nonunitary transformation, the Schrodinger equation reads 

H 

' A - i8 2V* 
2V -A + i5 

where 5 = (j a — 7&)/2, j a , lb being decay rates for upper and lower levels respectively; 

A = luq — v, ojq = (E a — Eh) , V = n ■ £, and /x is the electric dipole moment. To compare 

the geometric phase (IST)]) with that found in [3] we set x = 9ft V(t), y = $sV(t), z = A/2 

and e = 5/2. Then the geometric phase fl66l) being written as 

/ A — i5 \ , , 

7 = tt 1 _ . 68 

V ^/\2V \ 2 + (A -1S) 2 J 

coincides with the result obtained by Garrison and Wright [3] ■ 




)l 


U ) 




\ U2 / 



(67) 
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Remark. Comparing (1941) with Eqs. (lo"Tj) and (I60p . we conclude that both the Dirac 
complex monopole and hyperbolic monopole can be realized in the four-dimensional 
parameter space x, y,z,e G M 4 . The brief classification of the monopole structure is 
given in the Table 1. 



Monopole classification 


Dirac monopole 


£ = 


Complex monopole 


£ ^ 


Complex Dirac monopole 


£ 7^ 0, £ = const 


One-sheeted hyperbolic monopole 


£ ^ o, z = 0, x 2 + y 2 - £ 2 > 


Two-sheeted hyperbolic monopole 


£ ^ o, z = 0, x 2 + y 2 - £ 2 < 



Table 1. Monopole structure of the two- level dissipative system. 



4. Geometric phase and quantum evolution in vicinity of diabolic and 
exceptional points 

Since the adiabatic approach cannot be applied in the neighborhood of degeneracy, here 
we consider non-adiabatic generalization of the complex Berry phase. Let \u(t)) and 
(u(t)\ be solutions of the Schrodinger equation and its adjoint equation: 

i^\u(t)) = H\u(t)), (69) 

-i^(u(t)\ = (Z(t)\H (70) 

where we assume as usual the normalization condition (u(t)\u(t)) = 1. For an arbitrary 
evolution of the non-hermitian quantum system the complex geometric phase 7 = 7t — 7d 
is given by Eq. fl30l . and we have 

i , (u(t)\u(0)) . /•* , 
7= - In ),; ' \ t( +i / (u(t)\u(t))dt. (71) 



2 (u(0)Kt)> J 
This result can be adopted to calculate the geometric phase over the complex Bloch 
sphere as follows. Let n(t) = (sin a cos ft, sin a sin/?, cos a) be a unit complex Bloch 
vector defined as n(t) = (u(t)\er\u(t)) . One can observe that the Bloch vector satisfies 
the following equation: 

dn/dt = n x n, n(t) = Tr (H(t)tr), (72) 

and, as shown in the Appendix A, the complex geometric phase can be written as 

7 = — / (1 — cos a) ft dt + arctan — — ' -J^.. - — — — — . (73) 

1 2j y 1H cot(a//2) cot(ai/2) + cos{ft f - fa) V ; 

The integration is performed along the unique curve n(i) on the unit sphere S 2 , joining 
the initial point n(0) = n^ = (sinctj cos fti, sin«j shifts, cosctj) and the final point 
n(r) = Wf = (sin «j cos/3/, sinct/ sin/3/, cosa/). 
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Under a cyclic quantum evolution with the period T the Bloch vector describes a 
closed curve C on the complex two-dimensional sphere S 2 , and we have n(t + T) = n(t). 
The associated complex geometric phase, being half of the complex solid angle enclosed 
by C, 

l=-\ j{l - cosa)df3 = ~Q(C), (74) 

is known as the complex Aharonov- Anandan phase [U \57\ . 
Consider a generic non-Hermitian Hamiltonian 

H = ^l + ln.<r, (75) 

where 1 denotes the identity operator. Let \m) be a given initial state, then the 
solution of the Schrodinger equation (1751) can be written as \u(t)) = U(t)\ui) and 
(u{t) = {u^U^it), where 

U{t) = ( cos — - i sin y Cl ■ tr) e^ 2 (76) 
U-\t) = ( cos — + i sin — ft • ^V^ 2 (77) 



2 2 

Cl being the complex unit vector and fi = (fi • fi) 1 / 2 . 

Let \ui) and \v,f) = \u(t)) be initial and final states, respectively. Denoting 
the associated adjoint states by (ui\ and (uf\, we compute the transition amplitude 
Wi) — > \ui) and \ui) — > \uf) as 

T u =(u l \ U (t)\u l )=( K cos^-zsm^n l -Ciy- iXot / 2 (78) 

?>< = (uf\U(t)\ui) = ( cos % cosy -isiny n/i ■ ^)e iAo< / 2 (79) 

where cos^/j = rij = (■Uj|cr|it i ) and n/j = (w/-|<r|«i). 

The computation of the time dependent Bloch vector yields 

n(t) = cos fit rij + cos x(l — cos fit)f2 + sin fit (fi x rij), (80) 

where \ is the angle between the vectors rij and fi, so that is cosx = rij • Cl. 

Of the special interest is the case, when 9?(f2-f2) =0. Denoting by fio = ( | $~2 - 1 ) 1 / 2 , 
we obtain 

n(t) = cosfiotni + cosxo(l - cosfi t)^- + sm J^ ot (fi x rij), if Cl ■ Cl > (81) 

fio fio 

Cl sinh fi t 

n(t) = cosh fi t rij - cosxo(l - coshfi t) — + — x n*), if Cl ■ Cl < (82) 

fio fio 

where cos^o = n « ■ Cl/Q . At the exceptional point given by fi = and Cl = Cl e , both 
regimes yield 

t 2 

n(t) = rii - t(ni x n e ) + -(rij • n e )O e . (83) 



Similar consideration of the transition amplitude yields 
Tu = ( cos 2* - i sin ^ • e - iA °'/ 2 
2>< = ( cos cos f - i sin f n fi ■ e iA °'/ 2 



if O ■ > (84) 
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and 

Tu = ( cosh f - i sinh e~ iXot / 2 ) 

V oJ v iffl-n<0 (85) 

T fi = ( cos 6 fi cosh f - % sinh f • §-)e tX ° t/2 j 

Thus, if fi ■ f2 > we obtain coherent evolution of the quantum-mechanical system, 
and, if f2 ■ Q < we have incoherent one. At the exceptional point both regimes yield 



T- 



t 



(86) 

T/i = (^cos% - n/i • ft e Je* A ° t/2 (87) 
The complex geometric phase can be derived from the formula (1711) . or equivalently, 





Figure 3. Left panel: real part of geometric phase vs. a = 5R(n; • fi e ) and i 
(5R(nj • r2 e ) = 1). Right panel: K7 vs. b = 3(n j; • and t (3?(n 4 • ft J) = 0). 



using (|73p . Employing (1711) we obtain the time dependent geometric phase as 

fit i 1 + i cos x tanf 
7(*) = — cosx+ -In- : frr. (88) 



In the vicinity of the exceptional point we have 
7(*) 



t i 1 + ill; • il e I , „ N 

m • n e - + - in : — . ^ 2 + o(fi 2 ). 



1 - in.i • Q e | 



(89) 



2 2 

It follows from here that the real part of the geometric phase K7 (t) behaves like a 
step-function at the point to = 2/|3(nj ■ f2 e )| (see Figs. EH]). In particular, we obtain 

0, K(n, • n e ) -»• 0, t -»• t - 
-tt/2, 9&(n* • Cl e ) -f +0, t -f t + (90) 
tt/2, 3?(iij ■ ft e ) -»• -0, t -> t + 



K 7 (t ) = lim lim 9? 7 (t) 
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Figure 4. The graphic of 3?7 as function of a = SR(iij • fi e ) and b = • fl e ) (t = 2). 





) 









(91) 



4-1. Two-level atom driven by periodic electromagnetic field 

As an illustrative example we consider a two-level dissipative system driven by periodic 
electromagnetic field E(t) = §l(£(t) exp(iut)). In the rotating wave approximation, after 
removing the explicit time dependence of the Hamiltonian and the average effect of the 
decay terms, the Schrodinger equation reads [U EE] 

Ui \ _ 1 ( -«A + A - id 
u 2 )~2~\ 2V 

where A = (-y a + 7b)/2 with j a , 7 b being decay rates for upper and lower levels, 
respectively, A = u — u, u = (E a — Eb), 5 = (j a — 7b)/2, V — \x ■ £ and is 
the electric dipole moment. 

The choice £{t) = £ exp(iut) yields V(t) = Voexp(iut), where Vq — fj, ■ £o, and 
we assume further Vo > 0. The solution of Eq. fl9T|) with this choice of £ is well known 
and can be written as 

\ u (t)) = C 1 (t)e- i( "- iA) ' /2 |u T ) + C 2 (t)e l ^ +lX)t/2 \ Ul ) (92) 

where and |u-f) = ( o \ u l) = ( i ) denote the up/down states, respectively. In 
addition, \C(t)) satisfies the Schrodinger equation 

i 9 ^- = H r \C) (93) 

written in the co-rotating reference frame, where the Hamiltonian of the system takes 
the form 

H =±(*- U -* 2Vo ) (94) 
r 2\ 2V -A + u + i5 J [ ' 

and we find 

C x (t)\_( cos(fit/2) -icosxsin(fit/2) -i sin x sin (fit /2) 

C 2 (t) j~\ -ismxsm(Qt/2) cos(fit/2) + 2cosxsin(fit/2) 
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(95) 



where cos x = (A 



i5)/n, fi = (p 2 + (a 



10 



iS) 



2U/2 



2Vn. 



Without loss of generality, we can further confine out attention to the case e (p — 0, 
that implies Vq = [Vol - Passing on to the Bloch vector n(t) = (u(t)\cr\u(t)) we obtain 

' cos cut —sin cut o\ / sin 2 x + cos 2 % cos fit — cosx sin fit ~ sin 2%(1 — cos fit) 

sin cut cos cut cos x si n ^ cos fit — sin x si n fi^ 

^0 1 J \ \ sin2x(l — cos fit) sin x sin fit cos 2 x + sin 2 x cos fit 



n(t) 



Finally, one can show that n(t) satisfies the following equation 

da/dt = n'(t) x n 
where O' = (p cos out, p sin cut, A — iS) . 



(97) 



4-1.1. Cyclic evolution. Let n(t + T) = n(t) be he Bloch vector yielding a cyclic 
evolution of system over the complex sphere S 2 with the period T = 2tt/uo. Starting 
with its definition n = (u(t)\er\u(t)) , where \u(t)) and (u(t)\ satisfy the Schrodinger 
equation (19~TT) and is its adjoint equation, respectively, we find that the solution 



(w+(t)| =cos- e i( ^- 4 x)t/2 (u T I + sin ^ e -*<«-n+iA)*/2 1 ( u | 



K(t)) = cos I e 



(98) 
(99) 



yields 



(sin x cos cut, sin x sin cut, cosx)- (100) 

Note that n + (t) can be obtained as the periodic solution of the Bloch equation with 
fl = fi(sinx, 0, cosx)- The other periodic solutionis is given by 



n_ = — (sin x cos cut, sin x sin cut, cosx), 
and the related solution of the Schrodinger equation reads: 

sin 



u_(t)) = -sin I e -^- iA )*/> T > + cos ^e^- Q+iX W 2 \ Ul ) 



<«-(*)! 



- sin ? e «("+n-iA)t/2^ T | + cos X g -<( w -n+<A)t/2|^| 



(101) 

(102) 
(103) 

(104) 

/ P 2 + (A-uj-i5) 2 ' 

In the adiabatic limit, |cu/(A — i8)\ <C 1, the complex Aharonov-Anandan phase 
7 = 7_ + 27r reduces to the complex Berry phase 7^ obtained in [I]: 

A -«5 



2 * " 2 

The geometric phase derived from (|74p is given by 

A — cu — i5 



7± 



-7r(lT 



7ad = 7T 1 



\/ (2Vq) 2 + (A - i5) 2 



(105) 



where 2V = p. 
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In what follows we consider the behavior of the geometric phase 7 near the critical 
points, starting with the diabolic point. In this case $57 = and we have 

A ~ UJ N (106) 



7 = 7TI 1 



UJ 



1 



(107) 



This yields 

0, for p = 0, A - uj > 
2vr, for p = 0, A - uj < 

It follows that the geometric phase behaves as the step-function near the diabolic point, 
and at the diabolic point one has the discontinuity of 7 with the gap of 2tt (Fig. [5]). 





' 0.5 

'0.75 P 




Figure 5. Left panel: real part of geometric phase, 3?7 vs p and z = A — u>, in the 

vicinity of the diabolic point given by p = z = (5 = 0). Right panel: vs p and z 
nearby the exceptional point defined by p — 5 and z = (5 = 0.25). 

Once return to (11041) . we find that near the exceptional point, and for A = uj, real 
part of the geometric phase is given by 

, 7r, if p > 8 

** = { -(•■ , s < 108 > 



where the upper/lower sign corresponds to A — uj — ► ±0. Similar consideration of 
imaginary part of the geometric phase yields 

jr 0, if p< S 

9 H ,itp>6 < 109 > 

As can be observed in Fig. [51 the geometric phase has an infinite gap at the exceptional 
point. This is due to the coalescence of eigenvectors at the exceptional point. 
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4-1.2. Non-cyclic evolution. Let us consider the case in which the initial state is 
\u(0)) = | «-]•), which corresponds to the north pole of the Bloch sphere S 12 , and, hence, 



n(0) = (0,0,1). For non-cyclic evolution and the initial state chosen as n 
the explicit form of the time dependent solution of Eg. (1721) is given by 

' sin x cos x(l — cos fit) cos cut + sin \ srn &t sin out ' 
n(t) = sin x cos x(l — cos fit) sin ut — sin \ sin fit cos ut 
y cos 2 x + sin 2 x cos fit y 

The geometric phase derived from (173p is given by 



7 



fit tusin 2 ^-^ ^ N z' 1 + icosxtan^ 

k " Vl sin fit) + - In — 
; 2 1 



cosx - 



-(fit 



i cos x tan ^ 



(0,0,1) 



(110) 



(111) 



2 ""~' v 2fi 

where fi = (p 2 + (A — to — i^) 2 ) 1 / 2 . As depicted in Figs. Eland El the geometric phase 
has a singular behavior in vicinity of the exceptional point. At the latter, the imaginary 
part of the geometric phase shows the divergence and its real part behaves as step-like 
function. 




Figure 6. Left panel: imaginary part of the geometric phase, $57 as function 
of p and z — A — u> (t — 2/5, 8 = 0.5, ui = 1). Right panel: Q7 vs. p 
(z = 0,t = 2/6, 6 = 0.5, uj = 1). The divergence of S7 can be observed at the 
exceptional point {p = 5 = 0.5, z = Q). 



In the vicinity of the degeneracy point defined by fi = we have 
Zt . P V i 1+i f( 1 + i(fft 
7 2 6 + 2 l-if(l + l(f)*) +UlS1) ' lU2) 

where Z = A — uj — iS, and we set k — 1. If t 7^ 2/5, it follows that at the exceptional 
point, defined by Z = —i5 and p = 5, the geometric phase is given by 

tu5 2 t 3 5t i, 1 + f , . 

7 = z— + -n 1- 113 

' 6 2 2 1 - # v ; 
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0.49995 



Figure 7. Left panel: real part of geometric phase, 3?7 vs p and z = A — u> 
(t = 2/5 = 0.5, S = 0.5, u) = 1). Right panel: the same graphic nearby the exceptional 
point defined by p — 5 and z = 0. 



The case of t = 2/5 requires more careful analysis. Assuming t = 2/5 and inserting 
p = 5 into Eq. flll.2l) we obtain 

, 7r Auj A — uj j, e 2 |A — u;| ,„„ 
^ = ± 4-W + — -2 ll ' 2^ + i (A-L) +0( "> (U4) 
where the upper/lower sign corresponds to p — 5 — > ±0. At the exceptional point we 
have (Fig. El left panel) 

5ft 7 = ±T - ^, (115) 
In a similar way, assuming A — uj = 0, we obtain 

where the upper/lower sign corresponds to A — uj — > ±0. At the exceptional point we 
have (Fig. [SX right panel) 



4-1.3. Quantum evolution in vicinity of degeneracy. To study the tunneling process 
near degeneracy, we assume \u(t)) be a solution of Eq. (1911) with the initial state at 
t — chosen as |ttf), and the final state of the system at a later time t be or 
Then following [41. j, we compute the probability Pff (P^) that the system is in the state 
\ u t) (\ u l))i respectively, as 

P TT = |cos(fit/2) -icosx sin(fit/2)|V A * (118) 
P iT = |sin X sin(^/2)| 2 e- A * (119) 
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z p 

Figure 8. Left panel: graphic of vs z — A — to depicted for t = 2/8 and p = 8 
(8 = 0.5, uj = 1). At the exceptional point, defined by p = 8 and z = 0, the jump 
discontinuity is given by AK7 = ±7r/4. Right panel: vs. p (t = 2/8, z = 0. 
8 = 0.5, uj = 1). At the exceptional point = ±7r/2. 



In what follows we restrict ourselves by the case uj — A. Then, in according to the 
classification of the Tabl.l, the fictitious hyperbolic monopole emerges in the parameter 
space M 3 defined by the parameters of the system 9?Vo, QV , 5 G M 3 . 

There are two different regimes dependent on the relation between p and 5. For 
p > 5 we have one-sheeted hyperbolic monopole and coherent tunneling process 

p n = e - xt ( cos - — sin ) , (120) 



2 ft 2 

AT = e^sin 2 ^, (121) 

where f2o = |p 2 — 5 2 | 1 ^ 2 denotes the Rabi frequency. 

On the other hand, for p < 5, there is incoherent tunneling 
-\tf, Mot S . , fi n t\ 2 



P n = e -A* ^ cosh - — sinh -JL) , (122) 

P u = e--^sinh 2 ^ (123) 

and the associated monopole is the two-sheeted hyperbolic monopole. At the exceptional 
point we have Qq = 0, and both regimes yield 

2) 6 ^ ^ = ( 2 ) 
The Rabi oscillations is manifested as the quantum oscillation between up and down 

states and can be characterized by the following function: P(t) = P^ — P^ [59]. The 

computation yields 



H = 1 1 " "~ ) e~ Xt and P^=(-) e~ M (124) 



P(t) = e~ xt ( cos(fV) - — sin(fi t)) ,ifp>6 (125) 
P(t) = e ~ A ' C cosh(rV) - — sinh(fi t) ) , if p < 5 (126) 
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and at the exceptional point we obtain 

P(t) = e-*(l-!). (127) 

The Rabi oscillation function P(t) is plotted in Fig. [9] In addition, it is simple to show 
that in the absence of dissipation P(t) = cos(Oni). 





10 20 30 40 50 60 70 80 



Figure 9. Left panel: the Rabi oscillation P(t) as function oft for coherent tunneling 
( p > 6, Qo = 0.025, 6 = 0.1, A = 0.125). It is manifested as the quantum oscillation 
between the up and down states, and |it|). Right panel: P (t) versus t for 

incoherent tunneling ( p < S, CIq = 2, S = 0.1, A = 0.125). 



Once return to the geometric phase defined by Eq. flllll) . we obtain 



- - (hoi not - n t) -^ + ^ 00-6^ P ' 



20jj 



and 

7- ^3^ (sinhO t-M 



.ft 



1 fi + 5tanh^ 

2 n fi n -5tanh%^' 



if p < 5. 



:i28i 



(129) 



'0 - — 2 

It follows from Eq. (11281) that the real part of the geometric phase ffiry^t) has the 
jump discontinuity AK7 = =F7r/2 at the points 



2 / Q \ 
t n = — I irn ± arctan — ) , n = 0, 1, . . . 
IIq \ / 

with pulse duration At = t n +i — t n given by 



2tt , 



2 O 
— arctan — 

7T 



(130) 



(131) 



For the incoherent tunneling defined by Eq. (11291) the jump discontinuity AK7 = —n/2 
occurs at the point 

2 , n /On 



to = o; tanh " (f 



(132) 
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Figure 10. The real part of the geometric phase 3?7(r) versus time r = 2wt/oj. Left 
panel (one-sheeted hyperbolic monopole): coherent tunneling (f2 = 2, S — 0.5, uj = 1). 
Right panel (two-sheeted hyperbolic monopole): incoherent tunneling (fio = 0.25, 
S = 0.5, w = 1). 




Figure 11. Hyperbolic monopole. The real part of the geometric phase 3?7(r) versus 
time r = 2irt/uj and p is depicted. Left panel: dissipative system S ( 5 = 0.5, 
uj = 1). Right panel: 3?7(r) is plotted in absence of dissipation (6 = 0, u> = 1). As can 
be seen from the plot, the pulses presented at the left are disappeared. 

At the exceptional point Q = we have t = 2/5 and the pulse duration At — * oo. 
In absence of dissipation (5 = 0) the pulses are disappeared. Indeed, we have 
to = 2/5 — ► oo and the pulse duration At — >• while 5t — > 0. The real part of the 
non-adiabatic geometric phase 3?7(r) as function of the time r = 2-nt/io is plotted in 
Fig. [10], and in FigHJJthe real part of the geometric phase ^R.j(r,p) versus r and p is 
depicted. 
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To conclude, we note that the transition emerged at the exceptional p = 5 between 
two tunneling regimes is the topological phase transition in the parameter space, which 
can be described as follows: two-sheeted hyperbolic monopole (p < 5) <-> one-sheeted 
hyperbolic monopole (p > 5). 

5. Conclusion 

In the present paper we considered the geometric phase and quantum tunneling in 
vicinity of diabolic and exceptional points. It has been shown the complex geometric 
phase associated with the degeneracy points is defined by the flux of complex fictitious 
'magnetic' monopole. In weak-coupling limit the leading contribution to the real part of 
geometric phase is given by the flux of the Dirac monopole plus quadrupole term, and 
the expansion for its imaginary part starts with the dipolelike field. Recently similar 
result has been obtained for a two-level spin-half system in a slowly varying magnetic 
field and weakly coupled to a dissipative environment [60] . 

We found that the real part of the complex geometric phase has a discontinuity 
jump at the exceptional point. We also have shown that the exceptional point is the 
critical point of the quantum-mechanical system, and at the latter the topological phase 
transition in the parameter space occurs. 

Studying tunneling process near and at exceptional point we found two different 
regimes: coherent and incoherent. The coherent tunneling is characterized by the 
Rabbi oscillations, also known as quantum echoes, and it has been shown that the 
dissipation brings into existence of pulses in the real part of the geometric phase. At 
exceptional point both tunneling regimes yield the quadratic dependence on time, that 
is in accordance with the results obtained in [HU [62] for some specific non-Hermitian 
systems. The decay behavior predicted by Eqs. (11201) - (11241) has been recently observed 
in the experiment with a dissipative microwave billiard [62] . 

Emerging of pulses in the geometric phase is a novel quantum phenomenon, which 
reflects the monopole structure of the system. This complex-valued phase effect may 
be detected by using the quantal dissipative interferometry [BUI 161] . Note, that such a 
strong coupling effect of the environment should take place in generic dissipative systems, 
since in the neighborhood of the exceptional point only terms related to the invariant 
subspace formed by the two-dimensional Jordan block make substantial contributions 
and therefore the iV-dimensional problem becomes effectively two-dimensional [18| 150]. 
We conclude by remarking that obtained results would be important in implementation 
of fault-tolerant quantum computation, and it is necessary to understand better the 
relation between geometric phase and decoherence to perform computation [65, 66J. 
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Appendix A. Geometric phase for general evolution on the complex Bloch 
sphere 

In this Appendix we derive from the general definition of the geometric phase ( l30l ) 
written for two-level system as 

the formula (!73|) for computation of the geometric phase in terms of the complex Bloch 
vector. Eq. (1A.3I) generalizes to a non-hermitian Hamiltonian the formula obtained 
by Zang and Wang for computation of nonadiabatic noncyclic geometric phase for the 
Hermitian two- level system |67j . 

Theorem. The complex geometric phase defined in Eq. (1A.1I) is given on the 
complex Bloch sphere S% by 

1 f T n 1 h 2 -n 2 n 1 ( sux(/3/ - A) \ , ao x 

jit) = — / at + arctan — — J — — — , (A. 2 

1K ' 2j l + n 3 Vcot(a / /2)cot(a i /2) + cos(/3 / -A); 

where integration is performed along the unique curve n(t) on the unit sphere S%, 

joining the initial point n(0) = (sin a, cos/?j, sin sin/3j, cosaj) and the final point 

n(r) = (sina/ cos/3/, sin a/ sin/3/, cosa/). 

Proof. In general form, for a two-level system in terms of column and row vectors we 
have 

!«(*)>= f 6(J) V (u(t)\ = mXt))- (A.3) 

After some algebra and using the definition of the Bloch vector n(t) = (fy(t)\cr\ty(t)) , 
we find 

n\(t) = ab + ab, ^(t) = i{ab — ab), n^{t) = aa — bb (A. 4) 

From here, setting n(t) = (sin a cos /3, sin a sin/?, cos a), we obtain 

~ . a a _ ifj o a / \ r~ \ 

ao = sin — cos — e , aa = cos — , (A. 5) 

2 2 2 v 1 

~i ■ a a iR ,7 ■ 2 a i a ^\ 

ab — sin — cos — c , bb = sin — . A.o 

2 2' 2 v ; 

Next, denoting by |uj) and |ii/) initial and final states, respectively, we can write 
the total phase as follows: 

rt 2 V<«(0)l«(*)>/ 2 V<fiil«/>/ 2 ^2 i a / + 6 t 6 / l 
Then, applying (|A.5[) and (IA.6I) . we obtain 



_ » / cot(a//2)cotK/2) +e^ \ i /o.a/. 
2 Vcot(a//2)cot(ai/2) + e Wf-Pi)) 2 V 

a/Si/ 
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This yields 

/ sm(3f — 0i) \ i ( T ( da da\ , . 

7t = arctan - — / \ . P 7 — — — + - / . A.9 

Vcot(a / /2)cot(a i /2)+cos(/3 / -A)/ 2Jo \a a J ( J 

Since the dynamical phase 

o " r Jo 

we obtain 

sin(/3/ - Pi) 



ld =- l l (u(t)Au(t))dt = -i [ (dd + bb)dt, (A.10) 



7 = It ~ Id = arctan 



cot(a//2) cot(«j/2) + cos(/5/ - fa) 



i I j (Sd + 66) + - f | - - ) ) r//. ( A. I I } 







Using the relations Eqs. (1A.4I) - (1A.6j) . we find 



(1 - cos a)/? = n ^ 2 = H2hl = -2, f (fid + 66) + i f ? " - ) ) ■ (A.12) 



l + n 3 \ 2\da 

Then, inserting this result into (jA.llj) . we obtain 



7(r) = -- f T Hlil2 ~ H2hl dt + arctan f ? - , ^LzM ) . (A.13) 

V ; 2 J Q l + n 3 Vcot(a / /2)cot(a J /2) + cos(/3 / -A)/ 

□ 

Corollary. The geometric phase can be calculated by the following integral 

7(t) = -- / (1 - cosa)/M + arctan f 7 sm(/? / - fa) — ^\ 

rW 2 7 1 ^ Vcot(a / /2)cot( Q;i /2) + cos(/3 / -/3i)y 1 ' 

References 

[1] M. V. Berry. Quantal phase factor accompanyng adiabatic changes. Proc. R. Soc. London A, 
392:45 - 57, 1984. 

[2] Frank Wilczek and A. Zee. Appearance of gauge structure in simple dynamical systems. Phys. 

Rev. Lett, 52(24):2111-2114, 1984. 
[3] Y. Aharonov and J. Anandan. Phase change during a cyclic quantum evolution. Phys. Rev. Lett., 

58(16):1593-1596, 1987. 

[4] J. C. Garrison and E. M. Wright. Complex geometrical phases for dissipative systems. Phys. Lett. 
A, 128:177 - 181, 1988. 

[5] V. F. Weisskopf and E. P. Wigner. Berechnung der natiirlichen Linienbreite auf Grund der 

Diracschen Lichttheorie. Z. Physics, 63:54 - 73, 1930. 
[6] V. F. Weisskopf and E. P. Wigner. Uber die naturhche Linienbreite in der Strahlung des 

harmonischen Oszillators. Z. Physics, 65:18 - 29, 1930. 
[7] C. Jung, M. Miiller, and I. Rotter. Phase transitions in open quantum systems. Phys. Rev. E, 

60(1):114-131, 1999. 
[8] I. Rotter. Dynamics of quantum systems. Phys. Rev. E, 64(3):036213, 2001. 
[9] I Rotter. A continuum shell model for the open quantum mechanical nuclear system. Reports on 

Progress m Physics, 54(4):635-682, 1991. 
[10] M. V. Berry and M. R. Dennis. The optical singularities of birefringcnt dichroic chiral crystals. 

Proc. R. Soc. London A, 459:1261 - 1292, 2003. 
[11] M. V. Berry. Physics of nonhcrmitian degeneracies. Czech. J. Phys., 54:1039, 2004. 



Complex magnetic monopoles, geometric phases and quantum evolution 



27 



[12] M. V. Berry. Quantum adiabatic anholonomy. In U. M. Bregola, G. Marino, and G. Morandi, 
editors, Anomalies, Phases, Defects, page 125, Naples, 1990. Bibliopolis. 
M. V. Berry. Two-state quantum asymptotic. Ann. N. Y. Acad. Sci., 755:303 - 317, 1995. 
M V Berry. Proximity of degeneracies and chiral points. J. Phys. A: Math. Gen., 39:10013-10018, 
2006. 

Xiao-Chun Gao, Jing-Bo Xu, and Tic-Zheng Qian. Invariants and geometric phase for systems 

with non-Hermitian time-dependent Hamiltonians. Phys. Rev. A, 46(7):3626-3630, 1992. 
W. D. Heiss. Global and local aspects of exceptional points. Ann. Henri Poincare, 17:149, 2003. 
W. D. Heiss. Exceptional points of non-Hermitian operators. J. Phys. A: Math. Gen., 37:2455 - 
2464, 2004. 

W. D. Heiss and H. L. Harney. The chirality of exceptional points. Eur. Phys. J. D, 17:149, 2001. 
F. Keck, H. J. Korsh, and S. Mossman. Unfolding a diabolic point: a generalized crossing scenario. 
J. Phys. A: Math. Gen., 36:2125 - 2137, 2003. 

F. Keck and S. Mossman. Stark resonances for a double 6 quantum well: crossing scenarios, 
exceptional points and geometric phases. J. Phys. A: Math. Gen., 36:2139 - 2153, 2003. 

G. Dattoli, R. Mignani, and A. Torre. Geometrical phase in the cyclic evolution of non-Hermitian 
systems. J. Phys. A: Math. Gen., 23:5795-5806, 1990. 

Y. Aharonov, S. Massar, S. Popescu, J. Tollaksen, and L. Vaidman. Adiabatic measurements on 

mctastable systems. Phys. Rev. Lett., 77(6):983-987, 1996. 
S. Massar. Applications of the complex geometric phase for mctastable systems. Phys. Rev. A, 
54(6):4770-4774, 1996. 

A. Mondragon and E. Hernandez. Berry phase of resonant states. J. Phys. A: Math, and Gen., 
29:2567- 2568, 1996. 

G. Lindbland. On the generators of quantum dynamical semigroups. Commun. Math. Phys., 
48:119, 1976. 

D. Ellinas, S. M. Barnett, and M. A. Dupertuis. Berry's phase in optical resonance. Phys. Rev. 

A, 39(7):3228-3237, 1989. 
Dan Gamliel and Jack H. Freed. Berry's geometrical phases in ESR in the presence of a stochastic 

process. Phys. Rev. A, 39(7):3238-3255, 1989. 
A. Carollo, I. Fuentes-Guridi, M. Franca Santos, and V. Vcdral. Geometric phase in open system. 

Phys. Rev. Lett., 90(16):160402, 2003. 
M. B. Plenio and P. L. Knight. The quantum-jump approach to dissipative dynamics in quantum 

optics. Rev. Mod. Phys., 70(1):101-144, 1998. 
A. Carollo, I. Fuentes-Guridi, M. Franca Santos, and V. Vedral. Spin-1/2 geometric phase driven 

by decohering quantum field. Phys. Rev. Lett., 92(2):020402, 2004. 
Z. Fang et al. Anomalous Hall Effect and Magnetic Monopoles in Momentum-Space. Science, 
302:92-95, 2003. 

P. Zhang, Y. Li, and C. P. Sun. Induced Magnetic Monopole from Trapped A- Type Atom. Eur. 

Phys. J. D, 36:229 - 233, 2005. 
Rajendra Bhandari. Singularities of the mixed state phase. Phys. Rev. Lett., 89(26):268901, 2002. 
F. D. M. Haldane. Berry Curvature on the Fermi Surface: Anomalous Hall Effect as a Topological 

Fermi-Liquid Property Phys. Rev. Lett, 93:206602-1 - 206602-4, 2004. 
J. Frenkel and S. H. Pereira. Coordinate noncommutativity in strong non-uniform magnetic fields. 

Phys. Rev. D, 69:127702-1 - 127702-3, 2004. 
C. M. Savage and J. Ruostekoski. Dirac monopoles and dipolcs in ferromagnetic spinor Bose- 

Einstein condensates. Phys. Rev. A, 68:043604-1 - 043604-6, 2003. 
S. Murakami and N. Nagaosa. Berry phase in Magnetic Superconductors. Phys. Rev. Lett., 

90:057002-1 - 057002-4, 2003. 
Eytan Barouch and Barry M. McCoy. Statistical mechanics of the xy model, ii. spin-correlation 

functions. Phys. Rev. A, 3(2):786-804, Feb 1971. 
M. V. Berry and M. Wilkinson. Diabolical points in the spectra of triangles. Proc. R. Soc. London 



Complex magnetic monopoles, geometric phases and quantum evolution 



28 



A, 392:15 - 43, 1984. 

[40] Barry Simon. Holonomy, the Quantum Adiabatic Theorem, and Berry's Phase. Phys. Rev. Lett., 
51(24):2167-2170, 1983. 

[41] Howard C. Baker. Non-Hcrmitian quantum theory of multiphoton ionization. Phys. Rev. A, 
30(2):773-793, 1984. 

[42] W. D. Heiss. Exceptional points their universal occurrence and their physical significance. Czech. 
J. Phys., 54:1091, 2004. 

[43] O. Latinne, N. J. Kylstra, M. Dorr, J. Purvis, M. Terao-Dunseath, C. J. Joachain, P. G. Burke, 

and C. J. Noble. Laser-induced degeneracies involving autoionizing states in complex atoms. 

Phys. Rev. Lett., 74(l):46-49, 1995. 
[44] M. Philipp, P. von Brentano, G. Pascovici, and A. Richter. Frequency and width crossing of two 

interacting resonances in a microwave cavity. Phys. Rev. E, 62(2):1922-1926, 2000. 
[45] Markus K. Oberthaler, Roland Abfalterer, Stefan Bernet, Jorg Schmiedmayer, and Anton Zeilinger. 

Atom waves in crystals of light. Phys. Rev. Lett., 77(25) :4980-4983, 1996. 
[46] C. Dcmbowski,H.-D. Graf, H. L. Harney, A. Heine, W. D. Heiss, H. Rehfeld and A. Richter. 

Experimental observation of the topological structure of exceptional points. Phys. Rev. Lett., 

86:787- 790, 2001. 

[47] T Stehmann, W. D. Heiss, and F. G. Scholtz. Observation of exceptional points in electronic 

circuits. J. Phys. A: Math. Gen., 37(4):7813-7819, 2004. 
[48] V. I. Arnold. Geometric Methods in the Theory of Ordinary Differential Equations. Springer, 

New York, 1983. 

[49] P. M. Morse and H. Fcshbach. Methods of Theoretical Physics. McGraw-Hill, New York, 1953. 
[50] O. N. Kirillov, A. A. Mailybaev, and A. P. Seyranian. Unfolding of eigenvalue surfaces near a 

diabolic point due to a complex perturbation. J. Phys. A, 38:5531 - 5546, 2005. 
[51] A. A. Mailybaev, O. N. Kirillov, and A. P. Seyranian. Geometric phase around exceptional points. 

Phys. Rev A, 72:014104, 2005. 
[52] U. Guenther, I. Rotter, and B. F. Samsonov. Projective Hilbert space structures at exceptional 

points. J. Phys. A: Math. Theor., 40:8815-8833, 2007. 
[53] J. Anandan and L. Stodolsky. Some geometrical considerations of Berry's phase. Phys. Rev. D, 

35(8):2597-2600, 1987. 

[54] Peter Levay and Barnabas Apagyi. Algebraic scattering theory and the geometric phase. Phys. 

Rev. A, 47(2):823-830, 1993. 
[55] Roman Jackiw. Three elaboration on Berry's connection, curvature and phase, volume 5 of 

Advanced Serries in Mathematical Physics, chapter 1, pages 29-41. World Sci., Singapore, 

1989. 

[56] Luc Vinct. Invariant Berry connections. Phys. Rev. D, 37(8):2369-2372, 1988. 

[57] Shih-I Chu, Zhong-Chao Wu, and Eric Layton. Density matrix formulation of complex geometric 
quantum phases in dissipative systems. Chem. Phys. Lett., 157(1, 2):151-158, 1989. 

[58] Willis E. Lamb, Rainer R. Schlicher, and Marian O. Scully. Matter-field interaction in atomic 
physics and quantum optics. Phys. Rev. A, 36(6):2763-2772, 1987. 

[59] Chang-Qin Wu, Jian-Xin Li, and Dung-Hai Lee. Phase transition in a two-level-cavity system in 
an ohmic environment. arXiv.org:cond-mat/ '0703238, 2007. 

[60] R. S. Whitney, Y. Makhlin, A. Shnirman, and Y. Gefen. Geometric Nature of the Environment- 
Induced Berry Phase and Geometric Dephasing. Phys. Rev. Lett., 94:070407, 2005. 

[61] C. A. Stafford and B. R. Barrett. Simple model for decay of superdeformed nuclei. Phys. Rev. 
C, 60(5):051305(R), 1999. 

[62] B. Dietz, T. Friedrich, J. Metz, M. Miski-Oglu, A. Richter, F. Schafer, and C. A. Stafford. Rabi 
oscillations at exceptional points in microwave billiards. Phys. Rev. E, 75(2):027201, 2007. 

[63] Erik Sjoqvist, Aran K. Pati, Artur Ekert, Jeeva S. Anandan, Marie Ericsson, Daniel K. L. Oi, 
and Vlatko Vedral. Geometric phases for mixed states in interferometry. Phys. Rev. Lett., 
85(14):2845-2849, 2000. 



Complex magnetic monopoles, geometric phases and quantum evolution ... 29 

[64] Erik Sjoqvist. Quantal interferometry with dissipative internal motion. Phys. Rev. A, 
70(5):052109, Nov 2004. 

[65] V. Jones, J. A. and Vedral, A. Ekert, and G. Castagnoli. Geometric quantum computation using 

nuclear magnetic resonance. Nature, 403:869 - 871, 2000. 
[66] Shi-Liang Zhu and Z. D. Wang. Unconventional geometric quantum computation. Phys. Rev. 

Lett, 91(18):187902, 2003. 
[67] Shi-Liang Zhu, Z. D. Wang, and Yong-Dong Zhang. Nonadiabatic noncyclic geometric phase of a 

spin-i particle subject to an arbitrary magnetic field. Phys. Rev. B, 61(2):1142-1148, 2000. 



